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Abstract 

■ Assuming grand unified theory (GUT) and supersymmetry, we propose a simple 
model which can consistently accommodate the masses and mixings for quarks and 
leptons. The grand unified group is SO(10), and 10, 120, and 126 representations 
are introduced for the Higgs superfields which give masses to the quarks and leptons. 

\£) , The differences of masses and mixings between the quarks and the leptons are 

attributed to the Higgs boson structure. Below the GUT energy scale, the model 

■ is the same as the minimal supersymmetric standard model except its inclusion 
of dimension-5 operators for the small neutrino masses. The renormalization group 
equations of the independent parameters for the Higgs couplings with the quarks and 
leptons are given explicitly to connect the two energy scales of GUT and electroweak 
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theory. 
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1 Introduction 



Implications of non-vanishing neutrino masses are accumulating from the ex- 
periments for solar [1] and atmospheric [2] neutrinos which show neutrino os- 
cillations. This is the first experimental evidence that suggests physics beyond 
the standard model (SM). The possibility of oscillation has also been exam- 
ined by the neutrinos from nuclear reactors. The negative result of CHOOZ [3] 
gives constraints on the masses and mixing allowed for the neutrinos. The re- 
sult of KamLAND [4] confirms the oscillations of the solar neutrinos. It is 
a task for the extension of the SM to accommodate these experimental re- 
sults. In particular, the extreme smallness of the masses and the largeness of 
the generation mixing angles, both of which are observed in the experiments, 
should be accounted for naturally. 



Preprint submitted to Elsevier Science 



7 February 2008 



From a theoretical point of view, it is reasonable for the SM to be extended 
under grand unified theory (GUT). Furthermore, supersymmetry may have 
to be introduced in order that the SM be consistently embedded in GUT 
models. On the other hand, some of the GUT models, such as those based on 
SO(10) [5], automatically contain right-handed neutrinos. If they have large 
Majorana masses, the lightness of the neutrinos could be naturally understood. 
The plausibility of supersymmetric GUT models as the extension of the SM 
are then further strengthened. 

One tough obstacle confronts GUT models, concerning the generation mixing 
of the leptons. The lepton mixing is described by the Maki-Nakagawa-Sakata 
(MNS) matrix as the Cabibbo-Kobayashi-Maskawa (CKM) matrix does the 
quark mixing. Contrary to the CKM matrix, the observed neutrino oscillations 
show that the off-diagonal elements of the MNS matrix are not suppressed. 
However, since the leptons and the quarks are contained in the same rep- 
resentations of the grand unified group, similar magnitudes are deduced, by 
simple consideration, for corresponding elements of the two matrices. A suc- 
cessful GUT model is required to naturally implement the difference of mixing 
between the quarks and the leptons. Although there are various solutions ad- 
vocated for the problem [6], rather contrived schemes have been invoked. 

Aiming at a natural explanation for the generation mixings of the leptons and 
the quarks, we present a model base on SO(10) and supersymmetry in which 
the difference is simply attributed to the Higgs boson structure [7]. This model 
includes 10, 120, and 126 representations of the SO(10) group for the Higgs 
bosons which give masses to the quarks and leptons. The 120 representation 
becomes the origin of the different generation mixings. Large Majorana masses 
for the right-handed neutrinos are generated by a vacuum expectation value 
(VEV) of the 126 representation. Below the GUT energy scale, this model 
is described as the minimal supersymmetric standard model (MSSM) with 
dimension-5 operators composed of the SU(2) doublet superfields for the left- 
handed leptons and the Higgs bosons. The coefficients of the Higgs couplings 
with the quarks and leptons evolve between the energy scales of GUT and 
electroweak theory. We obtain the renormalization group equations for the 
eigenvalues of the coefficient matrices and the independent arguments of the 
MNS and CKM matrices. Quantitative analyses are then made on the neutrino 
oscillations. It is shown that the lepton mixing observed is compatible with 
the quark mixing in certain ranges of the model parameters. 

In SO(10) GUT models the Higgs bosons for the quark and lepton masses must 
be contained in 10, 120, or 126 representations. The direct product for one of 
these three representations with two spinor 16 representations for the quarks 
and leptons becomes an SO(10) singlet. The model with two Higgs fields of 
10 and 126 has been studied extensively [8]. However, all the experimental 
results for the neutrino oscillations are only marginally accommodated. The 
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10 and 126 representations make the Higgs couplings symmetric on the gener- 
ation indices, and their SU(2)-doublet components couple to both the quarks 
and the leptons. On the other hand, the 120 representation leads to antisym- 
metric Higgs couplings, and the SU(2)-doublet components which couple to 
the quarks may be different from those to the leptons. By including 120, the 
correlation of the masses and mixings between the quarks and the leptons 
becomes more flexible. 

In sect. 2 we summarize the interpretation of neutrino oscillations in the frame- 
work of the model with three generations for leptons. In sect. 3 the renormal- 
ization group equations for the MSSM with the dimension-5 operators are 
obtained explicitly in terms of independent parameters for the coefficients of 
Higgs couplings with quarks and leptons. In sect. 4 the GUT based on SO(10) 
and supersymmetry is discussed. After deriving in general the relations of 
Higgs couplings for the quarks and the leptons, a plausible model is speci- 
fied. In sect. 5 we make numerical analyses of the masses and mixing of the 
neutrinos. Conclusions are given in sect. 6. 



2 Neutrino Oscillations 

The deficit of neutrinos from the sun and the anomaly for neutrinos from 
the atmosphere can be explained by neutrino oscillations. The atmospheric 
neutrino observation gives that the mass-squared difference and the mixing 
angle between the /z-neutrino and a certain neutrino are given by 

Am 2 atm = (1 - 8) x l(T 3 eV 2 , sin 2 29 atm > 0.85. (1) 

There are several oscillation scenarios for the solar neutrino problem. However, 
a most likely solution is given by a large mixing angle between the e-neutrino 
and a certain neutrino under the Mikheyev-Smirnov-Wolfenstein effect, sug- 
gesting the ranges 

Am 2 sol = (2 - 10) x 10~ 5 eV 2 , sin 2 29 sol = 0.5 - 0.9. (2) 

The measurements by KamLAND are consistent with these results. On the 
other hand, by the experiment of CHOOZ the oscillation has not been observed 
for the e-neutrino, which excludes the combined region 

Am 2 chooz > 1 x lO^eV 2 , sin 2 29 chooz > 0.2. (3) 

These experimental measurements are taken into consideration in constructing 
our model. 
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We interpret the above results in the model with three generations of massive 
neutrinos. In this framework the strength of the W^-boson interaction for a 
neutrino and a charged lepton depends on the generations. The mass eigenstate 
v, L for the neutrino of the i-th generation is related to the eigenstate i>j of 
interaction with the charged lepton of the j-th generation by 

Vi = (V M Ns)iji>j, (4) 



where Vmns denotes the MNS matrix. The survival probability for the inter- 
action eigenstate i>i after run of distance L with energy E is given by 

(5) 
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J2 \( V MNs)ki\ 2 exp( 
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where m v \. represents the mass eigenvalue for the neutrino of the k-ih gener- 
ation. The mass-squared difference is hereafter written as 

Am?- = 77& - 771^-. (6) 



We assume that the neutrino masses are hierarchical, m v \ < m v2 < rn v3 , and 
not degenerated, similarly to the quarks and charged leptons. 

In the CHOOZ experiment, the phases in Eq. (5) are found to satisfy the 
relations (Am 21 /2E)L <C (Am 2 1 /2E)L ~ 1. The survival probability of the 
e-neutrino is expressed by 

P(u e - u e ) = 1 - 4|0W) 3 i| 2 (l - \(V M Nshi\ 2 ) sin 2 ~j^L. (7) 

Thus, the measured quantities are translated as 

Am 2 chooz = Am 2 .!, sin 2 26 chooz = 4|(lW) 3 i| 2 (l - IOW) 3 i| 2 ). (8) 



For the atmospheric neutrino oscillation, the phases in Eq. (5) satisfy the 
relations (Am 2 2l /2E)L < (Am\ 2 j2E)L and (Am 2 2l /2E)L < 1. The survival 
probability of the //-neutrino is given by 

A tyi 2 

Pfr - vy) = 1 - 4|(\W) 32 | 2 (1 - \(V MN s)s 2 \ 2 ) sin 2 -j^L. (9) 

The experimental results are expressed by 

Am 2 atm = Am 2 2 , sin 2 20 atm = 4|(\W) 32 | 2 (1 - \{V M nsU\ 2 )- (10) 



4 



For the solar neutrino oscillation, the relations |(Vm7vs)3i| 2 <C |(Vmats)2i| 2 
and \(V M Nshi\ 2 < |0W)n| 2 could hold from Eqs. (1), (3), (8), and (10). 
These constraints make it possible to evaluate the survival probability of the 
e-neutrino by 

P(u e - v e ) = 1 - 4|( W) 2 i| 2 (l - |( W) 2 i| 2 ) sin 2 -^L. (11) 

Therefore, we obtain the equations 

Aml ol = Ar/4, sin 2 26 sol = 4\(V MNS hi\ 2 (l ~ IOW) 2 i| 2 ). (12) 

These evaluations are used to discuss the predictions of the model. 



3 Energy evolution 



We assume that physics below the GUT energy scale is described by the 
MSSM with the dimension-5 operators composed of the superfields for left- 
handed leptons and Higgs bosons. This assumption is fulfilled in the model 
which is presented afterward. The superpotential relevant to the quark and 
lepton masses is given by 



w = w 1 + w 2 , 

W l = rflE x x Q l D cj + rfi H 2 x Q l U cj + 7#ifi x VE cj + H.c, 
W 2 = \^ ij H 2 x UH 2 x LP + H.c, 



(13) 
(14) 



where H\ and H 2 stand for the superfields for the Higgs bosons with hyper- 
charges — 1/2 and 1/2, respectively. Superfields are denoted by Q\ U a , and D cl 
for the quarks and U and E c% for the leptons, in a self-explanatory notation, 
with i being the generation index. The SU(3) group indices are understood. 
The dimension-5 superpotential W 2 yields the operators 



L = 



2 



V 



V 



+ H 2 ^ 




















U 









+ H.c, 



(15) 



where 0^- 2 and (p H2 represent the scalar components of the superfield H 2 . 
At the electroweak energy scale where the SU(2)-doublet Higgs bosons have 
non- vanishing VEVs, these dimension-5 operators become tiny Majorana mass 
terms for the left-handed neutrinos. The other operators arising from W 2 cause 
negligible effects. 
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The coefficient matrices rjd, r] u , r] e , and k are diagonalized by unitary matrices 
Uls and UrS. The CKM matrix for the quarks and the MNS matrix for the 
leptons are defined by 



V C KM = U?Ui, Uf Vu U^^^, Ui T Vd U d R=V^ (16) 
V M ns = U?UI U?kUI = k d , Uf Ve U e R=Ve, (17) 

where rjf, r/®, 77^, and k d stand for diagonalized matrices. Taking the VEVs 
of electroweak symmetry breaking for positive, the diagonal elements of rjf, 
r]®, and k d should be positive, while those of 77^ should be negative. The 
quarks and leptons then have positive masses. Those diagonal elements are 
expressed by i] di , r] ui , r] ei , and k { . 

A 3x3 unitary matrix has nine independent parameters. For the explicit ex- 
pression of Vckm or Vmns, we adopt the following parametrization [9]: 

V = P+Vo P_ j (18) 



^ S1S2C3 + e l<5 cic 2 cis 2 c 3 - e l5 sic 2 s 2 s 3 ^ 



V = 



S1C2C3 - e td cis 2 cic 2 c 3 + e td sis 2 c 2 s 3 



V 



-S1S3 -C1S3 c 3 

P+ = diag(exp(ipi), exp(ip 2 ), 1), 

P- = diag(exp(-ip 3 ), exp(-ip 4 ), exp(-ip 5 )), 



where q = cos#j and Sj = sin#j (i = 1,2,3). Without loss of generality, the 
angles Oi, 9 2 , and 6* 3 can be taken to lie in the first quadrant. At the electroweak 
energy scale, the numbers of the physical parameters for the CKM matrix and 
the MNS matrix are four and six, respectively. However, we need the general 
form of parametrization in Eq. (18) to discuss the energy dependencies of the 
CKM and MNS matrices. 

The values of 77^, rj®, 77^, k d , Vckm, and Vmns evolve depending on the 
energy scale. The renormalization group equations for these parameters and 
the gauge coupling constants of SU(3) x SU(2) xU(l) close on themselves at the 
one-loop level. We give the equations of the independent parameters, where 
the inequalities rf da < r]j 3 , r] 2 aa « 77^, and rf ea < 77^ (a = 1, 2) are taken into 
account. 

The gauge coupling constants: 

dgi _ o 93 dg 2 _ g 2 dg\ _ 33 g\ 



The diagonalized Higgs coupling coefficients: 
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9l ~ 3?& - - Ve 3 



13 
15' 



7 2 . 



Ids 



lu3 



(V M Ns) i3 \W e 3 



(20) 
(21) 

(22) 
(23) 



^ = -1^—3, (24) 

-f = -J^- 3, (25) 

^ = -^-3, (26) 

= ™ 

^^1 = ( a = 1-5), (28) 

where 8 q i and 5 q stand for the arguments in Vq for the CKM matrix, with 
Ci = cos 9 q i and Sj = sin 6^, and p ga represents an additional phase parameter 
in P + and P_. 
TTie MATS' mate 



dji 

d9n 
dji 
dOis 
dfj, 
dSi 
^ dfi 



// 



167T 2 



•S2C2C3 + C R )c 5 - (Bj - C^ss] , 



Ve3 



167T 2 



S2C2 



Ve3 



167T 2 



S3C3 



-A R S 2 + (Br + Cfl)cg 

-P^c 2 , + C R sl] , 



(29) 
(30) 
(31) 



16tt 2 



'■^i s 3 ~ (Br + Cr)ss— — 

S1C1 



S2C2C3 



-5/ ( c 5 — —s 2 c 2 c 3 + sjcj - C2 

V SlCl " , 
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p- 



dpjA 
dp 

dpi5 
dp 

A R - 

B R 

Br- 
C R - 



/ ™ / Ci S i 2 2 2 

+6/ c 5 s 2 c 2 c 3 - c 2 c 3 + s 2 

\ s i c i 

(A, - Bj)clsl + Cj(cl - slsl) 



(Aj - Ci)sls\ + Bj(4 - 4si) 



[B R + Cr)s & — S 2 C 2 
Si 

C l 2 \ f C l 
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Viz r 


dp 


16n 2 - 
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dp 


16tt 2 L 
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167T 2 3 [ 



Si 



+ 5/ c 5 — s 2 c 2 - c 2 c 3 - C/ s 2 c 2 + s 2 c 3 
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Ve3 

16n 2 
- Bj 

V 2 e3 

16tt 2 
2kiK 2 



(Sfl + Cr)s s — s 2 c 2 

Cl 



Si 



Cl 



C5 — s 2 c 2 + c 2 c 3 ) +Ci( c s — s 2 c 2 - s 2 c 3 



Si 



Si 



Cl 



s\ [Bjc 2 2 + C lS 2 2 



K 



r\ / \ ^1 ^2 

2 cos2(pi - p 2 ) + — ~~ 



K 



2 ' 
K 2 



2k±k 2 



K 



K 2 



sin2(pi - p 2 ), 



2k 2 k 3 . «1 + «I 

2 C0S 2 P2 + — " 
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3 
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2«3«i k| + k 2 
2 cos 2pi + — 
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2 ' 

K 3 — «i 



/l 'J 



sin2pi, 



(32) 
(33) 
(34) 

(35) 

(36) 
(37) 



where 9u and Si stand for the arguments in Vq for the MNS matrix, with 
q = cos 9n, Si = sin 9n, c$ = cos Si, and sg = sin Si, and p [a represents an 
additional phase parameter in P + and P_. The mutual dependencies of the 
independent parameters in energy evolution are manifestly seen, thanks to the 
explicit expressions in terms of the parameters themselves. For instance, the 
MNS matrix receives large quantum corrections if some of the neutrino mass 
coefficients «j are roughly degenerated [10]. 



4 Model 



The grand unified group of our model is SO(10). Its spinor 16 representation 
contains all the superfields for quarks and leptons of one generation, for both 
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left-handed and right-handed components. The right-handed neutrinos are 
therefore naturally incorporated. The decomposition of the direct product for 
two 16's is given by 16 x 16 = 10 + 120 + 126. The Higgs superfields which 
give masses to the quarks and leptons must be assigned to 10, 120, or 126 
representations. We introduce one superfield for each representation. 

The superpotential relevant to the quark and lepton masses are given, in the 
framework of SU(3)xSU(2)xU(l), by 



W = rf j 
+ e ij 



H[ Q x + UE cj ) + H[ Q x (Q l U cj + UN cj J 

H\ w x (Q l D cj + VE cj ) + ^#?fo x (Q iL>cj - 3^ l E cj ) 



+2H 5 120 x UN* + -j=Hfl x 

+ C j [H§g x (Q l D cj - 3VE cj ) + x (Q l U cj - 3VN cj ) 
+VQL iT H^U + ^H^N ci N cj ] + H.c, 



(38) 



where superfields N ct for the right-handed neutrinos appear in addition to 
the superfields Q\ U c \ D a , L\ and E a for the quarks and leptons. Higgs 
superfields are denoted by if's with upper and lower indices showing transfor- 
mation properties under SU(5) and SO(10), respectively: Hf , Hf 20 , Hf^o, and 
Hj^ are SU(2) doublets with hypercharge -1/2; Hf , Hf 20 , Hf 20 , and Hj^ 
are SU(2) doublets with hypercharge 1/2; Hj^r is an SU(2) triplet; and Hj^ 
is an SU(2)xU(l) singlet. Each superfield has been normalized. Owing to the 
antisymmetric property of 120, the coupling Hf 20 x Q l U cj does not appear. 
The coupling constants rfi and C u are symmetric for the generation indices, 
while e* J are antisymmetric. 

We can see from Eq. (38) the characteristic of the 120 representation. Any 
SU(2) doublet in 10 or 126 couples both the quark and the lepton superfields. 
As a result, in Eqs. (13) and (14), the Higgs coupling coefficients for the 
leptons become related to those for the quarks. On the other hand, for 120, 
four SU(2) doublets (VSHf 20 + ifj? )/2, (Hf 20 - VSHfl )/2, Hf 20 , and Hf 20 
couple to Q l D cj , L l E cj , L l N cj , and Q l U cj , respectively. The Higgs coupling 
coefficients could be less constrained. 



The SU(2)-doublet Higgs superfields for electroweak symmetry breaking are 
composed of the superfields in 10, 120, 126, and some other representations. 
Among the possible linear combinations of SU(2) doublets with the same 
hypercharge, only one doublet should be kept light to satisfy the unification 
of the gauge coupling constants for SU(3)xSU(2)xU(l). The two doublets 
with hypercharges —1/2 and 1/2 assume the role of the Higgs superfields Hi 
and H 2 in Eqs. (13) and (14). The other linear combinations must have large 
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masses and decouple from theory below the GUT energy scale. We express 
the Higgs superfields by 

H x = (Cj) u ff J, + (Cl) 12 Hf 20 + (Ct) 13 H? 20 + (Ct) 14 H§. + (39) 
H 2 = (CDn^o + (Cl) 12 H\ w + (Cl) w H« i + (C]) 14 ^ + (40) 

where C\ and C 2 represent unitary matrices. Some components of Hi and 
H 2 belong to the representations different from 10, 120, and 126, which are 
denoted by the ellipses. For instance, one superfield of 126 is included in 
the model. Its SU(5)-singlet scalar component has a large VEV to cancel the 
VEV of H^, and the VEVs of the auxiliary D fields for SO(10) are kept 
small. This 126 representation contains SU(2) doublets, which may become 
the components of Hi and H 2 . 

The matrices C\ and C 2 should be determined by the Higgs potential at the 
GUT energy scale. However, we take them as independent parameters in this 
paper, assuming that an appropriate Higgs potential could be constructed. 
Since the Higgs potential contains also the fields which are to break SO (10) 
correctly down to SU(3)xSU(2)xU(l), it is very complicated to analyze the 
whole potential. The Higgs potential is also supposed to induce the well-known 
split between the light SU(2) doublets and the heavy SU(3) triplets, as well as 
the split between the light SU(2) doublets, H x and H 2 , and the other heavy 
SU(2) doublets. 

The superpotential in Eq. (13) is now determined. The coefficient matrices are 
given by 



Vd = r](Ci)n + e[(Ci) 21 + ^(Ci) 3 i] + C(Ci) 4 i, (41) 

Vu = v(C 2 )n + ^<C 2 ) 3 i + C(C 2 ) 4 i, (42) 

Ve = v(Ci)u + e[(Ci) 2 i - V3(Ci) 3 i] - 3C(C 1 ) 4 i. (43) 

The dimension-5 superpotential in Eq. (14) is induced by the interactions of 
N ci in Eq. (38), 

W = 7]i j H 2 x VN cj + VQC, ij H^N ci N cj + H.c, (44) 

Vu = v(C 2 )n + 2e(C 2 ) 21 - 3C(C 2 ) 4 i. (45) 



If the scalar component of Hj^ has a large VEV, the right-handed neutrinos 
and sneutrinos receive large masses. The mass matrix of the neutrinos is given 
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by 

M VR = 2VSv s (, (46) 



with < >= v s /V2- The mass-squared matrix of the sneutrinos is ex- 

pressed as M UR Ml R . Exchanges of these particles lead to an effective superpo- 
tential given in Eq. (14), 

K = -r h ,(M VR )- 1 T£. (47) 



The coupling L lT Hj^D in Eq. (38) could give Majorana masses to the left- 
handed neutrinos, if the neutral scalar component of has a non- vanishing 
VEV. However, this VEV has to be as small as the neutrino masses measured 
in experiments. Then, an extreme fine-tuning of the Higgs potential would 
become inevitable. We therefore discard this possibility, assuming that 
is heavy enough not to develop a non-vanishing VEV. 

The coefficient matrices r/a, i] u , and rj e in Eq. (13) and k in Eq. (14) at the GUT 
energy scale are related to each other through Eqs. (41), (42), (43), and (47). 
For independent parameters we can take rju + Vui Vd + Vdi e > Ci> and C 2 . Then 
the symmetric parts of the coefficient matrices for the leptons are given by 



Ve + V I = 

T\ — r 4 
t 4rir 4 

Vu + Vv = 



3ri + r 4 ji. 4 . T . 

ivd + v d ) + - _ r iVu + Vu), 

' 1 '4 

[Vd + Vd) + {Vu + 



M, 



n 



ri - r 4 



V3v s 

(Ci) 4 l 

(C 2 )n 
(Ci)n' 



ri - r 4 



ri 



ri - r 4 



(Vd + Vd) ~ 



n - r 4 



iVu + Vu) 



(48) 
(49) 

(50) 



r 4 



2 Ml 



(Ci) 4 l" 



The mass matrix is symmetric. The antisymmetric parts of the coefficient 
matrices are given by 



Vd-Vd= 2 



^1)21 + ^(^)31 



^u-^ = ^(C 2 )3ie, 

r ?e -r ?e T = 2[(C 1 ) 21 - >/3(Ci) 3 i 
^-^ = 4(C 2 ) 21 e. 



(51) 

(52) 

(53) 
(54) 
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Depending on the structures of Eqs. (39) and (40), the 120 representation 
could contribute exclusively to any one of i] d , r) u , rj e , and r\ v . 

We now make an assumption that the (\^3Hf 20 +Hf2 )/2 and Hf^o components 
in the Higgs superfields H 1 and H 2 , respectively, can be neglected, taking the 
equations (61)21 + (l/v^) (61)31 = (62)31 = 0. Then, the coefficient matrices 
for the quarks become symmetric, Vd = Vd > Vu = Vu ■ Adopting a generation 
basis in which the coefficient matrix for the up-type quarks is diagonal, we 
obtain the equations 

Vd = VS KM ridV^ KM , r] u = rfc. (55) 



The matrices rj d and r\ u are determined by their eigenvalues and the CKM 
matrix. On the other hand, the coefficient matrices for the leptons are given 
by 



Ve 
Vu 



M, 



3rx + r 4 4 

-Vd H Vu + 4e, 



ri - r 4 ri - r 4 

4nr 4 n + 3r 4 _ 
-Vd H Vu + 2r 2 e, 



ri — r 4 



ri - r 4 



2^ 



(Ci) 4 l 

e=(Ci)2ie, r 2 = 



1 



ri — r 4 ri — r 4 
(62)21 
(Ci) 2 i" 



(56) 
(57) 

(58) 



The matrices rj e , r/ u , and M VR are expressed as linear combinations of 77^, 
and e. With being diagonal, the matrix i] d is roughly diagonal simultane- 
ously. However, for the matrix e, only off-diagonal elements have non- vanishing 
values. The contribution of 120 may make the off-diagonal elements of r) e 
and/or r\ v non-negligible, which could enhance the generation mixing for the 
leptons. 



The independent model parameters at the GUT energy scale are given by 
the diagonal matrices and 77^, the CKM matrix Vckm, the ratio n, r 2 , 
and r 4 , the antisymmetric matrix e, and the right-handed neutrino mass scale 
vs/ (Ci) 4 i. At the electroweak nergy scale, the eigenvalues of r] d , r) u , and r\ e are 
known experimentally, if the ratio tan (3 of the VEVs for Hi and H 2 is given. 
The CKM matrix has been measured. The quantities obtained experimentally 
for the neutrinos are the mass-squared differences and the mixing angles. These 
observed quantities have to be accommodated by suitable values of the model 
parameters at the GUT energy scale. 
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Fig. 1. The mixing parameters and the ratio of mass-squared differences for the 
neutrinos at the electroweak energy scale for the parameter set (A) : i) sin 2 29 a tm , 
ii) sm 2 26 sol , iii) sin 2 26 chooz , iv) Am 2 sol /Am 2 atm . 

5 Numerical analyses 

We show numerically that the observed neutrino oscillations can be described 
by our model. Since the number of the model parameters at the GUT energy 
scale is large, a systematical analysis in the whole parameter space is compli- 
cated. Instead, we demonstrate the viability of the model by presenting two 
numerical examples within the ranges of real values for ri, r 2 , r 4 , and e. 

In Figs. 1 and 2 the mixing parameters sin 2 2# aim , sin 2 2# so ;, sin 2 29 c f looz , and 
the ratio of mass-squared differences Am 2 ol /Am 2 tm at the electroweak energy 
scale are shown as functions of the model parameter r 2 . The values of tan/?, 
ri, r 4 , and e are listed in Table 1, with the sets (A) and (B) corresponding to 
Figs. 1 and 2, respectively. The other model parameters rj®, rj^, and Vckm at 
the GUT energy scale are tuned to give the quark and charged lepton masses 
and the CKM matrix at the electroweak energy scale shown in Table 2. The 
CKM matrix can be expressed by the standard parametrization [11] and the 
CP-violating phase of this parametrization is denoted by 5sp- These masses 
at the electroweak energy scale for the quarks and charged leptons, as well 
as the magnitudes of the CKM matrix elements, are consistent with the ex- 
periments [12]. The CP- violating phase 5sp lies in the range allowed by all 
the CP violation phenomena observed in the K°-K° and B°-B° systems [13]. 
The outcomes in Table 2 do not vary with r 2 . The value of vs/{Ci)n deter- 
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Fig. 2. The mixing parameters and the ratio of mass-squared differences for the 
neutrinos at the electroweak energy scale for the parameter set (B): i) sin 2 29 a t m , 
ii) sm 2 26 sol , iii) sin 2 26 chooz , iv) Am 2 sol /Aml tm . 

Table 1 

The parameter sets (A) and (B) for Figs. 1 and 2, respectively. 



(A) (B) 



tan/3 


20 


30 


n 


-3.9 


-1.9 


r 4 


-7.3 


-5.0 


ei2 


-4.0 x 10~ 4 


-9.0 x 10~ 4 


ei3 


3.6xl0~ 3 


4.9xl0~ 3 


£23 


9.2xl0~ 3 


7.6xl0~ 3 



mines the scale of k and does not affect the presented four observables for the 
neutrino oscillations. 

We can see from Figs. 1 and 2 that the atmospheric and solar neutrino os- 
cillations, under the constraints from the CHOOZ experiment, are realized 
simultaneously for certain parameter values. In Table 3 the resultant values 
of sm 2 29 atm , sin 2 29 soh sm 2 29 chooz , and Am 2 ol /Am 2 atm are explicitly given for 
the set (A) with r 2 = 12.5 and for the set (B) with r 2 = 7.5. If the mass scales 
of the right-handed neutrinos are put at vs/{Ci)±i = 1.0 x 10 15 GeV in the 
set (A) and at i>s/(Ci)4i — 6.0 x 10 14 GeV in the set (B), the mass-squared 
differences are given by Am 2 ol = 2.5 x 10~ 5 GeV 2 , Am 2 atm = 5.5 x 10~ 3 GeV 2 
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Table 2 

The masses of the quarks and charged leptons (in unit of GeV) and the CKM 
matrix (its elements and CP-violating phase in the standard parametrization) at 
the electroweak energy scale. 



m t 


1.8xl0 2 


\(Vckm)i2\ 


0.22 


m c 


6.8X10" 1 


\{Vckm)iz\ 


0.0036 


m u 


2xl0~ 3 


\{Vckm)2z\ 


0.041 


m b 


3.0 


Ssp 


1.3 


m s 


9.3xl(T 2 






m d 


5xl0~ 3 






m T 


1.8 






"V 


l.OxlO" 1 






m e 


5xl(T 4 







Table 3 

The mixing parameters and the ratio of mass-squared differences for the neutrinos 
at the electroweak energy scale with ri = 12.5 and r2 = 7.5 for the sets (A) and (B), 
respectively. 





(A) 


(B) 


sin 2 29 atm 


0.95 


0.94 


sin 2 29 so i 


0.59 


0.86 


sin 2 29 chooz 


0.16 


0.0026 



Am 2 sol /Am 2 atm 0.0045 0.055 

and by Am 2 sol = 8.7 x 10~ 5 GeV 2 , Am 2 im = 1.6 x 10~ 3 GeV 2 , respectively. 
Note that the mixing angle sin 2 29 c hooz is predicted to be around the present 
experimental bound for the case (A), while to be much smaller than it for the 
case (B). 



6 Conclusions 

We have discussed the masses and mixings of quarks and leptons within the 
framework of GUT coupled to supersymmetry. In GUT models the Higgs 
couplings for quarks and leptons are closely related to each other. The large 
generation mixing for the leptons, which is observed experimentally through 
the neutrino oscillations, cannot coexist trivially with the small generation 
mixing for the quarks. Some natural explanation for the difference of mixing 
between the quarks and the leptons is sought. 
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To solve the problem of generation mixing, we have proposed a model with 
SO(10) and supersymmetry. This model is a simple extension of the minimal 
SO(10) model, and includes a Higgs superfield of 120 representation, as well as 
two Higgs superfields of 10 and 126 representations. The energy evolution of 
the Higgs couplings with quarks and leptons have also been taken into account 
between the energy scales of GUT and electroweak theory. 

By the simple enlargement of the Higgs sector, the masses and mixings of the 
quarks and leptons are consistently accommodated without invoking much 
contrived schemes. The 120 representation makes the Higgs couplings different 
between the quarks and the leptons. The small neutrino masses are traced back 
to large masses for the right-handed neutrinos and sneutrinos generated by 
126. The model parameters are constrained to correctly give the quark and 
charged lepton masses and the CKM matrix. In spite of these constraints, 
the experimental results for the neutrino masses and the MNS matrix can be 
described well in certain regions of the parameter space. 
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